A unified theory for the cuprates, iron-based and similar superconducting systems: 
application for spin and charge excitations in the hole-doped cuprates 
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A unified theory for the cuprates and the iron-based superconductors is derived on the basis of 
common features in their electronic structures, including quasi-two-dimensionality, and the large-(7 
nature of the electron orbitals close to E-^ (smaller-C/ hybridized orbitals reside at bonding and 
antibonding states away from E-p). Consequently, low-energy excitations are described in terms 
of auxiliary particles, representing combinations of atomic-like electron configurations, rather than 
electron-like quasiparticles. The introduction of a Lagrange Bose field is necessary to enable the 
treatments of these auxiliary particles as bosons or fermions. The condensation of the bosons results 
in static or dynamical inhomogeneities, and consequently in a commensurate or an incommensurate 
resonance mode. The dynamics of the fermions determines the charge transport, and their strong 
coupling to the Lagrange-field bosons results in pairing and superconductivity. The calculated 
resonance mode in hole-doped cuprates agrees with the experimental results, and is shown to be 
correlated with the pairing gap on the Fermi arcs. 
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The recent discovery of high-temperature supercon- 
ductivity (SC) in iron-based compounds, including pnic- 
tides and chalcogenides Q (referred to below as 
FeSCs), provides an opportunity to test the validity of 
high-Tc theories in correlated-electron systems. Simi- 
larly to the cuprates [6|, the FeSCs are derived from 
an undoped "parent" compound which is generally mag- 
netically ordered [ri-loj at low temperatures and becomes 
SC under electron- or hole-doping. Also, both systems 
are characterized by a layered structure and a quasi-two- 
dimensional electronic structure 

A variety of normal-state properties including, e.g., the 
transport properties (i.e. resistivity, Hall coefRcient and 
thermoelectric power) of both the cuprates p7|-|2l| and 
the FeSCs [1, [13, are characterized by a remarkably 
similar anomalous behavior. Also, in both systems the 
suppression of SC by a high magnetic field results in a 
zero-ternperature insulator-to-metal transition upon dop- 
ing [23, [IHl • Even though the pairing symmetry is dif- 
ferent in the cuprates [26i - t28] and the FeSCs |29l - [33 |. a 
resonant spin excitation, characterized by wave vectors 
around those of the magnetic order in the parent com- 
pound, exists in the SC state of both systems 33 35]. 

The approximate tetrahedral arrangement of the pnic- 
togen/chalcogen atoms around the iron atoms in the 
FeSCs is typical of covalent bonding, and thus consid- 
erable hybridization is expected between orbitals cor- 
responding to the two atoms. This is confirmed in 
electronic-structure calculations however, such 

hybridization is found in antibonding and bonding states 
which lie at least ~ 1 eV away from the Fermi level 
{Ep ) , while the states at the close vicinity of are non- 
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bonding and of almost a pure Fe(3(i) nature |12Hla |. 

Consequently, the intrasite Coulomb and exchange in- 
tegrals, corresponding to Wannier functions of the hy- 
bridized orbitals of the entire bands which determine the 
Fermi-surface (FS), magnetic moments, etc., may be not 
large in the FeSCs [sBl, resulting in itineracy and largely 
reduced magnetic moments On the other hand, 

due to the dominantly Fe(3c?) nature of the states at E^, 
their intrasite integrals are rather large ^36| and a large- ?7 
approach should be applied to study the physical prop- 
erties {e.g. transport and SC) derived from low-energy 
excitations. 

This aspect of the electronic structure of the FeSCs 
is different from that of the cuprates, where an entire 
band around Ej„ is believed to correspond to the large- 
U physics [2^, Is^] , and an insulating state of large gaps 
and magnetic moments exists in the parent compounds. 
Low-energy carriers are present in the cuprates due to 
doping, and since such carriers in both the cuprates and 
the FeSCs correspond to the large-C/ regime, a unified 
theory could be worked out for both of them. This theory 
should be valid also for other quasi-two-dimensional SC 
systems which are close to a magnetic instability, and 
have large- 1/ electrons at the vicinity of E^ . 

At the basis of this theory stands the observation 
that SC exists at stoichiometrics where the dynamics 
of the low-energy carriers dominantly involves fluctua- 
tions between two adjacent occupation numbers (n) of 
atomic-like configurations (3d") around the copper or 
iron atoms. In the cuprates [28, Wf\ these are fluctua- 
tions between effective Cu(3c?^) (hybridized with 0(2p) 
orbitals) and Cu(3d^'') configurations for electron doping, 
and between effective Cu(3d^) and Cu(3c?^) (obtained 
through Zhang- Rice- type hybridization with 0(2p) or- 
bitals) configurations for hole doping. In the FeSCs these 
are fluctuations between Fe(3c?^) and Fe(3c?^) configu- 
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rations for electron doping, and between Fe(3(i^) and 
Fe(3d^) configurations for hole doping. 

Such dynamics of carriers could be treated through the 
auxiliary-particle approach [s^. A configuration corre- 
sponding to an occupation number n is denoted by a(ri), 
and a combined orbital-spin index of an atomic-like elec- 
tron by rj. For notation simplicity, let a{n— l{ri)) be the 
configuration obtained by removing an rj electron from 
a{n) 3 rj. The operator creates an auxiliary par- 

ticle representing the configuration a{n) at site i (a two- 
dimensional approximation is applied of points Ri within 
a planar lattice which could be defined to contain one Cu 
or Fe atom per unit cell [ill 113 ) ■ 

The creation operators of electrons of spin-orbitals 77 
at sites i can be expressed as: 

^ X/ X/ "iQ(n)"m(n-l(?)))- (1) 
« Q(n)977 

They satisfy anticommutation relations of independent 
fermion operators under the following conditions: (i) the 
consequence of the large- C/ approximation that only the 
contribution of two adjacent values of n could be con- 
sidered in rhs of Eq. ([1]) is valid; {ii) the auxiliary par- 
ticles created by , ^ are either bosons for even n and 
fermions for odd n, or fermions for even n and bosons for 
odd n; {in) the following constraint is satisfied in every 
site i: 

V Vflt a -1. (2) 

n a(n) 

As was discussed above, two occupation numbers (n) 
are considered, including no (corresponding to the par- 
ent compound), and either tiq + 1 (for electron doping) 
or no — 1 (for hole doping). Let us denote by a, /3 and 7 
the configurations corresponding to the occupation num- 
bers no -l- 1, no and no — 1, respectively. Their creation 
operators at site i are denoted by: 

Auxiliary-particles created by s^^ are chosen as bosons, 
and thus those created by e\^ and are fermions. 

The Hamiltonian "H, applied to study low-energy elec- 
tron excitations, is based on intrasite one- and two- 
particle terms, and intersite one-particle terms. It is 
expressed in terms of the auxiliary-particle operators 
through Eqs. (I1I3I) . A grand-canonical formalism is ap- 
plied by including in the Hamiltonian terms correspond- 
ing to the chemical potential /i, and to a field of Lagrange 
multipliers \i (A = (A^)) associated with the auxiliary- 
particles constraint [Eq. The values of Ai and ji 
should be determined to yield the correct charge and con- 
straint in every site. % could be, formally, expressed as 



(using constraint-preserving term): 

n = W +H''+H'' + AH, (4) 

= ^(4-A).t^s,^, 

ia. 

j^i a' PIS' 
if 

+ EE [KU^^-'^Mp'^,p^,hlY+h.c.]}, 

= EEE 

+ h.c], 

AH = -E(^'-^)[E4e.. + E4^./3 

i a p 

7 

The Ai — A Lagrange field represents an effective fiuctu- 
ating potential which prevents, through AH, constraint- 
violating fluctuations in the auxiliary-particle site occu- 
pation (thus enabling the treatment of atomic-like elec- 
tron configurations as bosons or fermions). The effect 
of such a fiuctuating potential on these configurations is 
analogous to the effect of vibrating atoms on electrons. 
Consequently, similarly to lattice dynamics, the quanti- 
zation of the Xi — A field yields bosons. 

In the cuprates one often applies a one-orbital model 
[28I. [37! ■ under which there is one a configuration, cor- 
responding to a complete Cu(3c?^") shell, one 7 config- 
uration corresponding to a Zhang-Rice singlet, and two 
P configurations corresponding to the spin states of the 
orbital a =t and a =\. (also presented here as a — ±1). 
The present auxiliary-particle method then becomes the 
"slave-fermion" method applied in previous works by the 
author 'i^, '4l| . The parameters appearing in Eq. Q are 
then simplified to the intrasite and transfer (hopping) 
integrals: 

(R) = Aw = t\l(^) = ^II(R) = m- (5) 

In the FeSCs one needs at least three Fe(3(i) orbitals 
[lM3 (of the xz, yz and — symmetries) to study 
the electrons at the vicinity of E-^ , and there are numer- 
ous a, /3 and 7 configurations. The parameters appearing 
in Eq. (0]) are then derived from intersite transfer, and in- 
trasite one-particle, Coulomb and Hund's-rule exchange 
integrals |4^ . 

Within the large-JJ approximation, applied in the 
derivation of "H, it could be approximated by omitting 
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FIG. 1: A typical lagron spectrum in hole-doped cuprates; 
the minima correspond to striped structures [i^ - lsH - 



either and the a term in AH, or and the 7 term 
in AH, and applying a second-order perturbation expan- 
sion in H'^'^. This results in corrections to hopping and 
intersite exchange terms which are expressed, within a 
one-orbital model for the cuprates, as: 

At(R,R')^-M^, J(R)^ffl^, (6) 



U 



U 



and Eq. Q is approximately replaced, for hole-doped 
stoichiometrics, by: 

■i <T 

+ [t(R, - Rj) + ^ At(R, - Rfc, Rfc - Rj) 



Considering values of i(R) up to third- nearest- 
neighbor R, and of t{R, R') and J(R) in Eq. © for 
nearest- neighbor R and R', yields an expression for H in 
terms of the parameters t, t' , t" and J. Values of these 
parameters for hole-doped cuprates have been obtained 
in first-principles calculations j43l - l46| . Explicit expres- 
sions for % (and its terms discussed further below), de- 
rived on the basis of Eq. ((T]) , will appear elsewhere [43| . 

The Lagrange field bosons are referred to as "lagrons" . 
They are soft at wave vectors corresponding to major 
fluctuations of spin and orbital densities. A typical la- 
gron spectrum in hole-doped cuprates is presented in 
Fig. [T] it has soft modes at the points: 



Qm = Q + (5q,„, for m = 1, 2, 3 or 4, 



(8) 



where Q = (7r/a)(x + y) is the wave vector of the anti- 
ferromagnetic (AF) order in the parent compounds, and 
5(\m = '^Sqx or ± 5qy are modulations around it, corre- 
sponding to striped structures [48l - [5l| . 



The Sj^-field bosons are referred to as "svivons" . Their 
Bose condensation is manifested, at low doping levels, in 
AF order, in the cuprates [i^, and in a structural distor- 
tion and magnetic order, characterized by a spin-density 
wave (SDW), in the FeSCs ^i^]. At higher doping levels 
the Bose condensation of svivons is manifested in static 
or dynamical inhomogeneities, based on modulations of 
the low-doping order. 

When svivons are Bose condensed, an s^^ field oper- 
ator can be expressed as a sum of its "condensed" part 
{i.e. the nonzero (s^^)) and fluctuating part s^^ — (Sj^). 
Thus, the expression of an electron creation operator in 
term of products of auxiliary-particle operators, through 
Eqs. (|1I3I) . includes terms where either or h^^ are 
multiplied by a condensed part of svivon operators, and 
terms where they are multiplied by their fluctuating part. 
A "quasi-electron" (QE) is defined as the fermion created 
by a normalized approximation to an electron creation 
operator, where only the terms in its expression which 
include condensed parts of svivon operators are main- 
tained. 

The QEs represent hypothetical approximate elec- 
trons which do not introduce fluctuations to the inho- 
mogeneities resulting from the Bose condensation of the 
svivon field. Since QE states are expanded as combina- 
tions of auxiliary-particle fermion states created by ei- 
ther the efo, or the /i^^ operators, these auxiliary-particle 
states form a basis to the QE states, and could be referred 
to as QEs as well. 

Thus, the problem of SC in strongly-interacting elec- 
tron systems is treated in terms of an auxiliary space 
consisting of three types of coupled "particles": (i) bo- 
son svivons which represent combinations of atomic-like 
electron configurations of the parent compounds, and 
their condensation results in static or dynamical inhomo- 
geneities; {ii) fermion QEs which represent combinations 
of such configurations with an excess of an electron or a 
hole over those of the parent compounds, and their dy- 
namics largely determines charge transport; {Hi) boson 
lagrons which represent an effective fluctuating poten- 
tial, enabling the treatment of the above configurations 
as bosons and fermions. 

Within this auxiliary space the pairing between the 
fermions through the exchange of bosons could be rigor- 
ously worked out in terms of coupled independent fields, 
in analogy to the electron and phonon fields within the 
BCS-Migdal-Eliashberg theory. The strong coupling be- 
tween QEs and lagrons, necessary for the constraint 
[Eq. ^] to be satisfied, results in high pairing temper- 
atures. If the same scenario were worked out as the 
pairing between electrons through the exchange of spin 
or charge fluctuations, generated by the same system of 
electrons, then two problems would have existed: (i) it 
is doubtful that such strongly-interacting electrons could 
be treated as quasiparticles; {ii) the coupled fermion and 
boson fields are not independent of each other. 

Svivon and QE spectra in hole-doped cuprates have 
been evaluated through a self-consistent second-order di- 
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FIG. 2: The absolute values of the svivon spectral functions below Tc for two typical hole-doped cuprates of different spin gaps; 
both the results for one condensate, and for their average over the four combined condensates are shown. 



agrammatic expansion [471 1 j where a mean- field treatment 
of y. in Eq. ([7]) is applied at the zeroth order. The ex- 
pansion is carried out on two Hamiltonian terms. One 
of them is AH which introduces svivon-lagron and QE- 
lagron coupling. Vertex corrections to it are negligible 
by a phase-space argument, as in Migdal's theorem, since 
the dominant contribution of the fiuctuating part of the 
constraint [Eq. ([2])] comes from a limited k-space range 
of the lagron spectrum around point Q (see Fig. [T]). The 
other term, "H', introducing QE-svivon coupling, is the 
contribution of the fluctuating part of the svivon op- 
erators to H. It is treated as a perturbation, and ap- 
proximated through a first-order expansion of the rhs of 
Eq. (O in terms like s^^ - (s-^) \4l\. 

Lagron spectra of the type presented in Fig. [T] deter- 
mine degenerate Bose-condensed svivon states, with en- 
ergy minima at points ±Qm/2, for one of the four values 
of m in Eq. Since there are four inequivalent values 
of Q/2 at ±(7r/2a)(x± y), the number of possible con- 
densates is eight. In the absence of symmetry-breaking 
long-range order, the system is generally in a combination 
of these states (refiecting fiuctuations between them). 
Tetragonal symmetry occurs when all the eight degen- 
erate states are combined, while orthorhombic symmetry 



breaking results in the combination of four of the eight 
states. The resulting stripe-like inhomogeneities p8l - [5l| 
(which resemble a checkerboard in the combination state) 
would be static or dynamical, depending on how close to 
zero are the spectrum minima. 

As they occur in Bose fields, the svivon spectral func- 
tions are positive at positive energies, and negative at 
negative ones. Their absolute values for two typical cases, 
of different nonzero spin gaps, in hole-doped cuprates 
below Tc (where the low-energy svivon linewidths are 
small) are presented in Fig. [21 Shown are the re- 
sults for the svivon condensate with energy minima at 
±[(7r/2a)(i + y) + ^Sqx], and the average of the results 
for the four condensates with minima at the vicinity of 
±(Tr/2a){x + y) (representing their combination), in ver- 
tical and diagonal directions around this point. 

The QE spectrum of hole-doped cuprates has been 
evaluated treating the fluctuations between the combined 
svivon condensates adiabatically, as is detailed in a sep- 
arate paper [52|. By the definition of electron creation 
operators in Eq. ([Ij, the electron Green's functions are 
obtained at the zeroth-order as sums of products of QE 
and svivon Green's functions. This results in the non- 
Fermi-liquid (non-FL) scenario of a distribution of con- 
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FIG. 3: The imaginary part of the spin susceptibility, corresponding to the svivon spectra for hole-doped cuprates below Tc 
presented in Fig. [21 the small- and large-spin-gap results demonstrate, respectively, the existence of an incommensurate [sj or 
a commensurate y3| resonance mode. 



voluted QE-svivon poles. It is shown [521 ^-^^^t multiple 
scattering of QE-svivon pairs introduces to the electron 
Green's functions additional FL-like electron poles, and 
thus the effect of both types of poles is reflected in various 
physical properties. 

The spin susceptibility (SS) of hole-doped cuprates has 
been evaluated, under an approximation where only the 
non-FL convoluted QE-svivon poles are considered [s^ . 
Linear-response theory has been applied on the basis 
of spin-flip processes, expressed by constraint-preserving 
terms of the form {sl^s^ _^s^j _^Sj^) , and thus deter- 
mined by the svivon spectrum. Results obtained for the 
imaginary part of the SS, in vertical and diagonal di- 
rections around k = Q, are presented in Fig. [3] They 
correspond to the two svivon spectra shown in Fig. [2l 
and since the svivon-system is in a combination state, 
the SS results are averaged over those of the four com- 
bined condensates. These results reproduce those ob- 
served in neutron-scattering measurements in different 
hole-doped cuprates. The larger spin-gap results cor- 
respond to the "commensurate resonance mode (RM)" 
[331] . and the smaller spin-gap results correspond to the 



"incommensurate RM" [33 |. 

If the constraint is imposed in any two sites i and j, 
the following equation should be satisfied in these sites 
in hole-doped cuprates [see Eqs. (|2I3I) ]: 

(Tcr' 

The terms in the Ihs of Eq. ([9]) are formally similar to 
the above spin-flip term applied for the derivation of the 
SS results presented in Fig. [31 Thus, a susceptibility- 
like function, referred to as the "constraint susceptibility" 
(CS) could be derived on the basis of either the svivon 
spectrum, through the Ihs of Eq. ^ , or the QE spectrum, 
through the rhs of Eq. ©. The results obtained for the 
CS on the basis of both spectra should agree with each 
other in order for the constraint to be satisfied, and this 
condition is the basis for the determination of the lagron 
spectrum, and of their coupling to the svivons and QEs. 
The CS represents the response of auxiliary particles, and 
not of electrons. However, it reflects, under certain con- 
ditions, an approximation to the response of the system 
to charge fluctuations which could be measured, e.g., by 
Raman spectroscopy [1^, [s^ . 
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FIG. 4: The imaginary part of the constraint susceptibility corresponding to the svivon spectra for hole-doped cuprates below Tc 
presented in Fig. (5] the low-energy peaks around k = 0, approximately, correspond to the integrated energies of the resonance- 
mode peaks shown in Fig. [S] these peaks also, approximately, correspond through Eq. ((9]) to the SC gap over the Fermi arcs 
(see discussion in the text). 



Results obtained for the imaginary part of the CS, on 
the basis of the Ihs of Eq. in vertical and diagonal 
directions around k = 0, are presented in Fig. S) They 
correspond to the two svivon spectra shown in Fig. [5J 
and evaluated similarly to the SS results presented in 
in Fig. 131 The major feature observed in CS results is a 
low-energy peak around k = at energies which, approx- 
imately, correspond to the energies of the k-integrated 
low-energy features of the SS at the vicinity of k = Q 
(thus the incommensurate or commensurate RM). 

Since the same CS results, as those presented in Fig.|4l 
should be obtained also on the basis of the QE spectrum 
through the rhs of Eq. and since they correspond to 
the SC state, the observed peak at k = should represent 
some kind of an average value of the QE gap below Tc- 
As is explained elsewhere fs^l, this gap has two contri- 
butions; one of them originates from Brillouin zone (BZ) 
ranges around the antinodal points, where a narrow peak 
(of energy e = at T = 0), lying between two humps, 
splits due to pairing below T*; the other contribution to 
that gap opens below Tc on the Fermi arcs (FAs) around 
the line of nodes. 



In the SC state there are both "normal" and "anoma- 
lous" (pair-correlation) QE Green's functions, and their 
contributions to the QE expression for the CS have op- 
posite signs [13). These contributions cancel each other 
for "gap-edge states", where e = 0, E = Ve^ + A-^ = A, 
and thus the fraction of both the particle and the hole 
states within the Bogoliubov states is ^[1 ± e/E] = i. So 
the QE-spectrum contributions to the CS peak at k = 
come from states where e 7^ 0. 

Consequently [E^l, the averaged QE gap which deter- 
mines the k = CS peak is lowly weighted around the 
antinodal points, and represents a value somewhat larger 
than the averaged QE gap on the FAs. Since the averaged 
electron FA gap is also somewhat larger than the QE FA 
gap (due to convolution with svivon states) , one expects 
a correlation between the values of this gap and the k = 
CS peak, and as was discussed above (see Figs. [3] and H]), 
also with the averaged RM energy. The electron FA gap 
has been measured through, e.g., the B2g Raman mode, 
and its value has indeed been found to be correlated with 
the RM energy [26l.[27l[33j. A correlation between the en- 
ergies of the Aig Raman mode and the RM [5^ has been 
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found to be partial |54|. The observed correlation of the 
FA gap with ~ 5k^T^ j26l.[27j is explained elsewhere [53 |. 
The fact that the average RM energy is lower when it 
is incommensurate (see Fig. [3]) explains the observation 
that Tc is lower in cuprates with an incommensurate RM. 

Even though the electronic structure of low-energy 
states in the FeSCs is based on more orbitals than in the 
cuprates, important physical conclusions could be drawn 
from one system to the other due to the formally com- 
mon Hamiltonian applied for both of them. Within a 
two-dimensional approximation, the lagron spectrum of 
the FeSCs is expected to differ from that of the cuprates, 
presented in Fig. [l] by replacing the minima positions 
from satellite points around Q, to satellite points around 
the two possible SDW wave vectors in the parent com- 
pounds: Qi = {Tr/a)x and Q2 = {'n/a)y 0, S ( O'" 
Qi = {'K/2a){x + y) and Q2 = {■K/2a){x — y) [9|), o r 
points close to them. Similarly to the cuprates [48145 1| . 
stripe-like inhomogeneities characterized by modulations 
due to the differences between the satellite points and Qi 
or Q2, could exist also in the FeSCs. 

The svivon spectrum in the FeSCs is expected to have 
analogous features to those of the cuprates, presented 
in Fig. [21 resulting in a resonance mode in the vicinity 
of Qi and Q2, below Tc, as has been observed [S^. In 
a separate paper [55], it is explained that QE pairing 
requires a sign reversal of the order parameter upon a 
shift of Q in the BZ, in the cuprates, and of Qi or Q2 



in the FeSCs. Due to their different FSs, this results in 
pairing symmetry of an approximate d^2_y2 type in the 
cuprates, and of an s± type (thus with different signs on 
different FS pockets) in the FeSCs. Thus, it is predicted 
that there are no Fermi arcs in the FeSCs, and that their 
RM energy is correlated with an averaged value of the 
SC gap, as has been observed f30l-l32i [ssj . 

It could be concluded that high- Tc SC occurs in 
quasi-two-dimensional strongly-interacting electron sys- 
tems due to the fact that low-energy excitations in them 
are described in terms of auxiliary particles, represent- 
ing combinations of atomic-like electron configurations, 
rather than electron-like quasiparticles. A Lagrange Bose 
field which must be introduced to enable the treatments 
of these auxiliary particles as fermions or bosons, serves 
as the pairing glue between the fermions. 
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